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AN EVALUATION OF FOUR 


SINGLE ELEMENT AIRFOIL ANALYTIC METHODS# 

R. J. Freuler and G. M. Gregorek 
General Aviation Airfoil Design and Analysis Center 
The Ohio State University 


SUMMARY 


A comparison of four computer codes for the analysis of two-dimensional 
single element airfoil sections is presented for three classes of section geo- 
metries. Two of the computer codes utilize vortex singularities methods to ob- 
tain the potential flow solution. The other two codes solve the full inviscid 
potential flow equation using finite differencing techniques, allowing results 
to be obtained for transonic flow about an airfoil including weak shocks. Each 
program incorporates boundary layer routines for computing the boundary layer 
displacement thickness and boundary layer effects on aerodynamic coefficients. 

Computational results are given for a symmetrical section represented by 
an NACA 0012 profile, a conventional section illustrated by an NACA 65A413 pro- 
file, and a supercritical type section for General Aviation applications typi- 
fied by a NASA LS(l)-0413 section. Experimental results from The Ohio State 
University 15 cm (6 in.) by 56 cm (22 in.) Transonic Airfoil Tunnel are also 
given. The cases presented include operating conditions at subsonic, sub- 
critical, and near critical or supercritical Mach numbers. The four codes are 
compared and contrasted in the areas of method of approach, range of applica- 
bility, agreement among each other and with experiment, individual advantages 
and disadvantages, computer run times and memory requirements, and operational 
idiosyncrasies. 


INTRODUCTION 


The General Aviation Airfoil Design and Analysis Center (GA/ADAC) was esta- 
blished at The Aeronautical and Astronautical Research Laboratory (AARL), The 
Ohio State University, under contract to NASA Langley Research Center in June 
1976. GA/ADAC offers a comprehensive service to the general aviation community 
in the form of airfoil selection and design and analysis work as well as con- 
sultation in the areas of wind tunnel testing and flight testing work. An im- 
portant feature of GA/ADAC is the large library of computer codes which has 
been established and is maintained at AARL. This computer program library re- 

#This work has been supported in part by NASA TIngley Research Center ”” 

Contract NAS1-14406. 
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location nn+^ de variety of airfoil related computer codes collected in one 
areas ^of < f 3 ^ 60 ’ wlth more than thirty computer codes in the 

Sa?vs?s airf ° 1 11 T lyai8 desi ^> roulti-element airfoil 

fo™2ce ^lysIs^aSbirSrSr" Pr0Peller aerodj ™” 10 “ d aa ° aa « a p«- 

at eva J ua J io V f the foUr CQm P^er codes most frequently used 

is t 5 ( “ >W » of two-dimensional single element airfoil sections 

d in t JJ 8 P?P er for three classes of airfoil section geometries- 

oTJr l dZl °H ly Ca ! bered ' and a ^t -cambered. TheoreticS predictions 
fl rp P .!l^ distributions and aerodynamic coefficients for the three airfoils 

taken in the 15 cm n 56 cm TranaoSe AirtoU 


The symbols used herein are defined in an appendix. 


COMPUTER CODES 


/. i The computer codes used for the comparisons are designated as follow*? • 
tailerSI^’ by F, 1 BaU f; P / ° arabadl ' m ’ D. Korn, and fl^eaon ^d S: 

^firencea 3 (2) v°^ 1 T’ by L ' A - Carla °" md plained in 

references j and 4 , (3) Smetana, by F. Smetana, D. Summey, N Smith and R 

Carden and documented in references 5 and 6; and ( 4 ) Eppler, by R. Eppler aAd 

D. Somers and soon to be documented in a NASA Technical Note 7 It should be 

mentions that tho versions of these counter codes in use at GA/ADACare min- 

thi tl t S Up ' t0 " da I' e as P° ssltl «J yet, in some instances, these versions are not 

code L cS^V 8 ' 11 f. the COdes ’ ” lth tha of the Smetana 

code, is constantly being refined and improved by the respective program 

the^ef f e ct s°of ^^aminar^h Ca f son code is currently being modified to include 

fied to itL«tp , t bol f y layer ^ the ^P ler cod ® soon be modi- 
tied to iterate on the boundary layer displacement thickness. 

given\erJ- f tSe S m^. 1 , 0 ure f c- + h %"f th ? d °, f ?? proaoh for the four “O' 165 >» 

X r«S?!r literature cited (refs. 1-6) contains the detailed explanations 

tion to the f!m°i e ^ I* transonic code that employs a finite difference solu- 

fJil The bo^Lr nV i 1SCl 5? te ? tial flow equation for a conformally mapped air- 
a ary iayer displacement is added iteratively to the airfoil 

effects The°Carlsnn ® ValU ? te airfoil section performance including viscous 
ects. The Carlson code is also a transonic code and is similar to 

arabedian, but it uses a finite difference solution to the full inviscid pot- 
ential flow equation for an airfoil in a stretched Cartesian coordinate svs?L 

iub^i?i 0f | f °V con f on " ally ma PPed airfoil. The Smetana program is a stJi^tly 
B n^ Crl ^ 1C + 1 a C v de a method of vorticity distributed over an airfoil 

approximated by a closed polygon. An iterative approach to bo^da“ Xyer 
effects is included. The approach used in the Eppler code is similar to the 

fonZ^ ’ bUt / td i f f erS in that the vorticity is distributed oJL Air- 
foil shape approximated by curved panels and that there is no iteration on the 
boundary layer displacement thickness. on the 


que of me ! h ° dS 1 0f approach for the lnviscid flow and the techni- 

2. ?! handling the boundary layer for each computer code is given in table 1 

The iterative procedure referred to in table 1 is that of: (l) obtainlm L 

JoT, i +? id f 10W , SOlution for the ori & ln al airfoil; (2) obtaining a boundary layer 
soiution based on the inviscid flow solution; ( 3 ) modifying the airfo^LaS 

by adding the boundary layer displacement thickness to the airfoil- ( 4 ) obtain 
ing an inviscid flow solution for the modified airfoi?; 2 d, rUatiS 
steps 2 through 4 until convergence criteria are satisfied. 


AIRFOIL SECTIONS 


tive r2,n+2 a T S ° f ai f f0 . il f ection geometries are included in these compara- 

m2ted Mr! m!2 f f ? ° J a ftion is re P resen ted by the well docu- 

mented NACA 0012 profile (ref. 7 and 8) shown in figure 1. To provide a suffi- 

eent2 er °, ordinates for computational purposes, computer 
generated coordinates for the NACA 0012 (ref. 9) were used with tS aWM! 
section being defined by 47 ordinates for both i^^Jd ^r spaces M 

h % S b6en made in tUs t0 a ™luate ST peSo^Sce of 

silts 2 a2foil°o22at C °f S Tf th respect t0 sensitivit y of computational re- 
least 35 to b de ^ S f ty ’ a sufflclent nuraber of coordinates (i.e. at 

least 3 ? to 40) has been used to provide consistent, reliable results Ordinate 

density has been distributed such that there is a higher densi2 of point! 22 
centrated In regions of greater airfoil curvature. Ms ~ 

quirement for the two subcritical codes because of the nature of the distri- 
buted vorticity methods of flow solution. distri- 

foii ?fif t I) f0r T L C ° n 2 nt + i0na i Se ^J° n " e Hl us Crated by a NACA 65A413 air- 
the method of 3rdi " ates for this section, 63 in number, were obtained by 

NACA 6A 22?! ! +? nCe 10 ; . Some dra & Prediction comparisons for another 
NACA 6A-series section, a 64A010, are also included. 

w . Th ® third c } ass of airfoil geometries investigated is the aft-camber*»d 
Whitcomb supercritical type section. In this paper, the results forTde^L- 

the 6 NASA S LSn^ o??? i0 ? ape ® lficaU y for general aviation applications, 

the £! 3 2 0il (f £’ 3) ; are Presented. The ordinates used for 

iefeSnce'n ' ° &S the GA(W) ’ 2 airfoi1 ’ those listed ln 


airfoil data summary comparisons 


.. Fo f ! h a, tW0 subc ritical codes, Smetana and Eppler, a comparison of the 

? ir !? ilS 2 ! emS '« airfo11 data summary Plots is 0 r?n2res?. As 
in for the NACA 0012, the comparison with experiments reported 

!olid f i?2 whIi l8 +h Ui 2 g °° d 4 N0t ® that the wind tunnel test is shown as the 

veitio2 Thl Ll^t I 8 given by the symbols ' a reversal of usual con- 

ventions. The Eppler code shows a break-over in the Cl (lift coefficient) 


/ 

i 


Y 


. J”L 


versus alpha (angle of attack) plot produced by using the predicted separation 
point to define an "effective" angle of attack. The Smetana code merely iden- 
tifies a predicted separation point, hut makes no attempt to compensate for the 
effects of separation. The free transition option was specified in the Smetana 
code for this comparison and the following two airfoil data summaries. The 
Eppler code always uses natural transition, although it does allow a variable 
roughness option. 

In figure 5, computational results are shown for a NACA 65A413 airfoil 
compared to the wind tunnel results for a NACA 65]_-412 of reference 7. 
codes may be observed to predict the laminar drag bucket for this NACA 6-series 

section. 

The airfoil data summary comparisons for the NASA LS(l)-0413> shown in 
figure 6, point out the difficulty Eppler has with some airfoils with regard to 
the angle of zero lift. This appears to be related to the lack of an iterative 
boundary layer solution and is more noticeable with supercritical typ®, blunt 
trailing edged airfoil shapes. The wind tunnel results are those of McGhee, et 
al. (ref. 12). 


PRESSURE DISTRIBUTION COMPARISONS 


The detailed pressure distribution comparison cases which follow include 
operating conditions at subcritical and near critical or supercritical Mach 
numbers. The computational results are compared with experimental results from 
The Ohio State University 15 cm (6 in. ) by 56 cm (22 in. ) Transonic Airfoil 
Tunnel. The OSU 6 x 22 wind tunnel is a low- interference transonic facility 
for airfoil testing over the Mach number range of 0.30 to 1.07 and a Reynolds 
number range of 2 to 15 million based on 15.24 cm (6 in. ) model chord (refs. 

13 and 14). The angle of attack used in the computational results Is the 
effective angle of attack, ct e ff, obtained from the set angle of attack in the 
wind tunnel corrected for wall effects according to the empirically derived 
relation: 


a eff = “set “ 0,17 C L 


It should be noted that for all the supercritical pressure distribution 
comparisons, the drag coefficient listed for Carlson has been omitted. Total 
drag as predicted by Carlson requires a wave drag correction to be applied that 
was not available to the authors at this writing. 


Figure 7 shows the comparisons for all the codes for the NACA 0012 airfoil 
at a Mach number of about 0.35 over a range of angles o f attack. The pressure 
distributions are in good agreement, though the Eppler code predicts a somewhat 
higher suction peak than the other codes at the higher angles of attack. Drag 
comparisons for the subcritical codes using the free transition option are 
quite good. The less accurate drag predictions by the two transonic codes is a 
result of attempting to simulate free transition by fixing turbulent boundary 
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layer transition a few percent chord in front of the Smetana predicted natur- 
ally occurring transition location. 

- C^P^isons for the NACA 0012 at an angle of attack of zero over a range 
of Mach numbers is given in figure 8. Both the Smetana and Eppler codes show a 

Mach 6 n2w Pr6 !}ji lower Pressures in the nose region at the near critical 
Mach number condition, and of course can not correctly predict the distribution 
at supercritical conditions. Both transonic codes identify the strength and 
location of the shock quite well. ^ 3110 

m „ Comparisons of the transonic codes and wind tunnel results for the NACA 

noilno?w U n e ftn rit J Cal c ?? dltlons ^ g iven in figure 9. The Mach number is 

? ?' 8 ° re fnlts are given for three angles of attack. Note that the 
arlson code appears to predict somewhat higher values of lift and a corres- 

Prediction of a shock located further aft on the airfoil. This appears 
b ® caused by an uncertainty in angle of attack in the Carlson code, with a 
trend toward results being obtained at a slightly higher angle of attack than 
^ J?? 11 Sngl ? 0f attack for many airfoils. Thus, in general. Carlson results 
Sb °? ld b ® ex f lined as Pressure versus lift coefficient, moment versus lift S- 

£ “ * tC -> taad of of attack. Since direct comparisons with wind 

tunnel angle of. attack were desired for this study, matching angle of attack 
as more convenient, so this approach has been used. Both codes indicate a 

1° re fJ er /°^ e Pressure on the aft upper surface than is observed in 
Q test !. for tMs airfoil (and most other airfoils as well). In 

thfTthP Jpd i Cr u Pa !! Cy ShOCk location be explained by the fact 

that the predicted local Mach number m front of the shock is in excess of 

wind tunnel ° k nUmb6r that poses diff iculty for both the theory and the 

For the transonic codes, careful selection of input parameters relating to 
ducI 6 2v nCe £ e d axatl °? factors are required to encourage the codes to pro- 
urt results for an airfoil when the free stream Mach number and 

lift coefficient exceed certain values. The empirical relationship below, 
suggested by Dr. R. Whitcomb, appears to describe these limiting values: 

M + t/c + 0.1 Cj, _> 0.92 

5 T er ? the free stream Mach number, t/c is the airfoil thickness ratio and 
OL is the lift coefficient. 

In figures 10 and 11, results are presented for the NACA 65A413 airfoil 

t S +J Cr J tiC v 1 J f ld slightly supercritical Mach number. This compari- 
son shows botn the Garabedian and Carlson codes over-predicting the lift 
Carlson's over-prediction could be related to the angle of attack uncertainty 
previously discussed but no consistent reason can be presented for Garabedian' s 
results, particularly for the generally higher pressures predicted on the lower 
surface (fig. 10). This characteristic in Garabedian occurs infrequently and 
may be circumvented by using the matching-lift-coefficient option (which 
essentially compensates for any uncertainty in angle of attack in either the 
wind tunnel or computer code). Also, the theoretical predictions of drag, 
though consistent with each other, are lower than the wind tunnel results. 
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(The wave drag contribution in figure 11 is less than 0.0002, aid thus the 
Carlson skin friction drag as shown approximates the total drag). 

Comparisons for the NASA LS( 1 )-0413 airfoil section are given in figures 
12 and 13. Figure 12 shows the computational and experimental results obtained 
for a nearly zero angle of attack through a range of Mach numbers from 0.45 to 
0.80. For this airfoil, agreement with the wind tunnel results at subcritical 
conditions is excellent for both the pressure distributions and the aerodynamic 
coefficients. It is interesting to note in figure 12c that the pressure 
distributions and shock locations predicted by Garabedian and Carlson differ 
noticeably, presumably for reasons mentioned earlier; yet in figure 12d the 
pressure distributions predicted by the codes are nearly identical. (The wind 
tunnel results in 12d may be influenced by the strong shock present at the 
condition illustrated). This comparison points out the uncertainties in angle 
of attack are dependent not only on the input airfoil but also on the specific 
input conditions as well. 

The results for the NASA LS(l)-0413 at a nominal Mach number of 0.72 over 
a range of angles of attack are presented in figure 13. In figure 13a, both 
transonic codes exhibit some interesting characteristics. Carlson, although 
the lift nearly matches the wind tunnel results, has difficulty properly de- 
fining the nose region on the lower surface. This may be due to the fact the 
Cartesian grid used does not place a large number of computational points near 
the leading (and trailing) edge. Although the Garabedian result accurately 
describes the lower surface nose region including the shock location, the lift 
prediction is too low as a result of over-predicting the pressure recovery on 
the upper surface. In figure 13b these same trends may be observed to a lesser 
degree. In 13c the results are more characteristic of the codes: Garabedian 

showing a reasonable lift and drag prediction with a slightly higher than wind 
tunnel observed pressure recovery over the trailing edge region; and Carlson 
•exhibiting results at an apparently higher effective angle of attack for the 
input angle of attack which matches the wind tunnel and the Garabedian results. 

In figure 14, comparisons are shown for Carlson with wind tunnel results 
(fig* 13c) by both matching angle of attack and selecting angle of attack which 
is matching the wind tunnel lift coefficient. Note that the expected excellent 
agreement of the pressure distribution and the aerodynamic coefficients with 
the wind tunnel test is obtained when lift coefficients arc matched. 


DRAG PREDICTIONS 


Drag coefficient predictions by the Garabedian code are generally consis- 
tent and accurate enough to enable use of the code to predict the drag rise 
characteristics and the drag divergence Mach number for most airfoils. The 
version of th? Garabedian cede in use at GA/ADAC for the past 18 months employs 
the latest wave drag calculation techniques and the fast Poisson solver for the 
subsonic region of flow which improves the rate of convergence (ref. 2). 
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figure 15 shows good agreement between the Garabedian code predicted and 
the OSU 6 by 22 wind tunnel observed total drag for Mach numbers well into the 
drag rise. It should be noted that the Reynolds number was not to be held con- 
stant in these results, but varied from 3.8 to 5.9 million. The angle of 
attack was nominally zero. Transition was specified as fixed at 0.075c for the 
Garabedian code which seems consistent for the LS(l)-0413, a turbulent flow 
airfoil by design, at these conditions. The sudden decrease in skin friction 
drag shown for this airfoil in the drag rise region is a result o? shock-in- 
duced separation due to the strong shock present. The relatively larger re- 
gion of computer predicted separated flow behind the shock on the upper sur- 
face of the airfoil results in a lower skin friction drag coefficient. 


Drag rise characteristics for the symmetric, NACA 64A010 airfoil section 
are given in figure 16. Excellent agreement between theory and wind tunnel is 
again observed. The boundary layer transition was fixed at 0.05c in both the 
wind tunnel and computer code. The Reynolds number varied from 3.5 million at 

Mach 0.5 to about 5 million at Mach 0.85. The angle of attack was held at 
zero. 


COMPUTER REQUIREMENTS 


All four airfoil analysis codes are run on the GA/ADAC computer facility 

AARL ’ The com Puter system is a dual processor system using Harris 
SLASH 6 and SLASH 5 processors. The SLASH 6, which is used for all of 
GA/ADAC s airfoil work, is a medium-sized, 24 bit word computer system with 
b4K words (192K bytes) of main memory. For purposes of comparison, the SLASH 
o is about an average factor of 8 times slower in heavy floating point FORTRAN 
programs than an IBM System 370 Model 168. Of special interest is the fact 
the calculations on the SLASH 6 are performed with over 11 decimal digits of 
accuracy while single precision on IBM mainframes affords approximately 7 deci- 
mal digits accuracy. The 11+ digit accuracy of the computer used at GA/ADAC 
is well suited for most scientific calculations including airfoil analysis, 
thus avoiding the necessity to maintain 15-16 digit accuracy like that of a 
GDC mainframe or double precision on an IBM machine. 

2 lists several computer related characteristics of the four codes. 
Of most Interest are the memory requirements and the run times for the programs. 
All the codes have been folded into the SLASH 6 such that the largest program 
requires 48K words . All programs are overlaid to varying degrees to reduce 
memory requirements. Per case run times, where a "case" is a calculation at 
one Mach number, one Reynolds number, and one angle of attack, are expressed in 
a normalized form. The single case run time used for normalization is that of 
the Smetana code. For the SLASH 6, time T is on the order of 90-100 seconds. 

A range of times is shown for the two transonic codes since convergence to a 
solution varies depending on whether the case is subcritical or supercritical 
and on a user-supplied convergence tolerance. The Carlson code has the long- 
est running time per case when results are carried to the fine grid, which was 
used in all cases for the previous compa. isons. The medium grid result of 
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time . ^ Smetan^code^a^no Snvercenee 0 ^^ 71 ^ reduction in computation 

gram-fixed number of inviscid flow/bo^d^layef i+eratioS ° n a pro ‘ 

verged solution and resulting in jTv ^i 0118 to achieve a eon- 

Eppler code is the most rapid o? /T , tlme per case ‘ ^ 

percentage of time T per case All nf i^ cd ® s ®? d typically requires a small 
plot capability. m of the codes ^ve some form of hard copy 


OBSERVATIONS AND CONCLUSIONS 


fo „ _ ° n , the r 68 ^ 8 presented here and the extensive exercise of these 

our single~e lenient airfoil analysis codes +he f*niirt , nr? < nrf _ . • 
marks can be made. ys cooes, the following observations and re- 

SMETANA CODE: is a subcritical code; uses vorticity distributed 

BwSSSSa^^ war - 

a curved STairfor ® Sub ° r jf tlcal code »* uses vorticity distributed around 

wide v g °° b ?? ndar3r layer rou tines, applicable over a 

wide range of Reynolds numbers; exhibits good performance at low Revnolds 

numbers; has a design mode, although the mode is difficult to use at first* 

effeots estimate ’ «*** 

. ° oefficient versus angle of attack; executes very quickly, resulting in 

KS Pe Bu? a n- C S U n 1 "? t C °M i “ d ' 0ntalns ™«VirlkL TX 

difficulty in £^£*5 th e 

18 related to n ° *•—£ 

THE GARABEDIAN CODE: is a transonic code; employs a finite different® 

airfoil f? th ! fUl1 ? nviscid Potential flow equation for a conformally mapped 
airfoil, iterates on boundary layer disnlappmpn+ ^ 

tiStTZXsT* “ “ “^S^S’-SE E2 SoT,T 

P de f reasonable wave drag estimates and car. be used for drag rise 
predictions; and has flexible input options, allowing to specify eithe- angle 

$j£&L°ZSi: , I i * imat , qe Un • But lt: has no ^ inar boundary layer & or^ e 
f t . e . r a ’ employs a boundary layer smoothing process which can tend 

si ° w do - s — « - 
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THE CARLSON CODE: is a transonic code; employs a finite difference solu- 

tion to the full inviscid potential flow equation for an airfoil in a stretched 
Cartesian coordinate system; iterates on boundary layer thickness; gives good 
pressure distributions; has easy to use design mode; and incorporates a massive 
separation prediction technique. But it: has no laminar boundary or transition 
criteria (but one is currently being added); does not have an input option for 
matching lift coefficient; exhibits an uncertainty in angle of attack; and 
needs improvement in prediction of the wave drag coefficient. 


APPENDIX 


SYMBOLS 

Measurements and calculations were made in the U.S. Customary Units. 
They are presented herein in the International System of Units (SI) with* the 
equivalent values given parenthetically in the U.S. Customary Units. 

a angle of attack, deg 

c chord 

Cjj drag coefficient 

Q 

L lift coefficient 

(^1 pitching-moment coefficient 

C p pressure coefficient 

K = 1024 

M Mach number 

RE Reynolds number 

T computer solution time per case 

"k/ c airfoil thickness-to-ehord ratio 
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TABLE 1 . COMPARISON OP POUR SINGLE ELEMENT AIRFOIL ANALYSIS CODES 
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Figure 4.- NACA 0012 airfoil-section characteristics. Computational results 
are for a Reynolds number of 6 million at a Mach number of 0.20. 
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Figure 5.- NACA 65A413 airfoil-section characteristics. Computational results 
are for a Reynolds number of 6 million at a Mach number of 0.20. Comparison 
is made to a 65.-412 airfoil of reference 7. 
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Figure 6.- NASA LS(1)-0413 airfoil-section characteristics. Computational 
results are for a Reynolds number of 6 million at a Mach number of 0.20. 
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Figure 7.- Concluded. 


LIFT 

□ .830 

0 .839 

A .907 
+ .891 

X .865 

MOMENT 

□ -.006 

0 .029 

A .010 
+ -.001 
X -.010 

DRAG 

a .0100 

0 .0085 
A .0090 
+ .0126 
X .0103 




LIFT 


□ 

.006 

o 

.001 

A 

0.000 

+ 

0.000 

X 

0.000 

MOMENT 

□ 

-.001 

0 

-.000 

A 

0.000 

+ 

0.000 

X 

0.000 


DRAG 

□ 

.0081 

0 

• 0047 

A 

.0051 

+ 

.0061 

X 

.0065 


(a) M = 0.575; RE = 4.68 million. 
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(b) M » 0.725; RE = 5.34 million. 
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(a) M - 0.808; RE = 6.12 million; ct = 0°. 
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Figure 9.- Comparison of computer-code predictions with wind-tunnel results 
for an NACA 0012 airfoil section at a supercritical Mach number over a 
range of angles of attack. 
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Figure 10.- Comparison of computer-code predictions with wind-tunnel results 
for an NACA 65A413 airfoil section at M * 0.517, RE = 7.08 million, and 
a = 2.90°. 



Figure 11.- Comparison of computer-code predictions with wind-tunnel results 
for an NACA 65A413 airfoil section at M ■ 0.700, RE ■ 8.20 million, and 
a = -0.06°. 
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(a) M - 0.454; RE - 3.75 million; a - -0.07 
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(b) M - 0.578; RE - 4.98 million; a - -0.07 . 

Figure 12.- Comparison of computer-code predictions with wind-tunnel results 
for an NASA LS(1)-0413 airfoil section at an angle of attack of zero over 
a range of Mach numbers. 
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Figure 14.- Comparison of results from the Carlson code obtained by attempting 
to match lift coefficient to wind-tunnel result. NASA LS(1)-0413 airfoil; 

M = 0.722; RE = 4.69 million; a = -0.09°. 
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Figure 15.- Comparison of the Garabedian code prediction and wind-tunnel result 
for the drag-divergence characteristics of the NASA LS(1)-0413 airfoil 
section at zero angle of attack. 
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Figure 16.- Comparison of the Garabedian code prediction and wind-tunnel result 
for the drag-divergence characteristics of the NACA 64A010 airfoil section 
at zero angle of attack. 



